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Abstract 

This paper reports a new phenomenological analysis of the neutrino burst detected from SN 
1987 A, and it reveals the presence of two mass eigenstates. The heavier mass eigenstate has 
mu = 21.4 it 1.2eV/(? , while the lighter one has = 4.0 it 0.5eV/(P . It is not the first paper to 
make such a claim, but it expands on a 1988 conditional analysis by Cowsik, and it attempts to make 
the evidence more robust through an improved statistical analysis, and through providing reasons 
why alternative explanations are unlikely. It also shows how the result can be made consistent with 
existing smaller electron neutrino mass limits from tritium beta decay and cosmology through the 
existence of a third tachyonic (superluminal) mass eigenstate. 

PACS numbers: 13.15+g, 14.60Pq, 14.60St, 14.60Lni 



I. INTRODUCTION 



The neutrinos from SN 1987A have been the subject of hundreds of papers, both theo- 
retical and phenomenological. jl| Some of these papers analyze the data to infer an upper 
limit on the neutrino mass, which ranges typically from 12 to 16 eV/c^.^, 3|. It is expected 
however that neutrinos from a supernova should include all three active neutrino mass eigen- 
states, therefore such a limit must represent some sort of average over whichever eigenstates 
are represented in the neutrino events recorded. There have in fact been at least three 
analyses of the SN 1987A neutrino data that have given evidence for a bifurcation into two 
groupings. I^-Ig I . The analysis by Loredo was based on Baysian statistics, jsl and it gave no 
physical interpretation to the two classes of events. The Roos paper based the bifurcation on 
neutrino mass, but he attributed the different mass values to the events recorded in the two 
detectors, rather than to two mass eigenstates. Finally, the Cowsik paper did suggest a bi- 
furcation associated with two mass eigenstates whose masses he found based on the relation 
between neutrino arrival times and their measured energies. Cowsik reported these eigen- 
states as 24 ± TeV/c^ and 4 ± leV/c^. We here extend his analysis to include the possibility 
of observed regularity arising in part from chance, or variations in their emission times at 
the source, since his conclusion was conditional on there being near instantaneous emissions. 
Additionally, we show how the result can be made consistent with existing electron neutrino 
mass limits, and other experiments. As with most SN 1987A analyses, we include neutrino 
events recorded in the Kamioka, and 1MB detectors, but exclude the 5 events seen in the 
Mont Blanc detector since they occured 5 hours before those detected by the other three 
detectors. Additionally, given its size and sensitivity that detector should have only seen 1 
neutrino from the supernova not 5. We also exclude 5 events seen in the Baksan detector 
because there are substantial reasons to believe that perhaps 2 or 3 of them may have been 
due to background. 

As noted in the paper on the Baksan neutrinos from SN 1987A.[3] 

• The 5 event burst they saw over a 10 s period was only slightly above their background, 
since there have been past instances when the Baksan detector recorded 3 "random" 
events in 10 seconds 

• The timing of the Baksan burst occurred 30 s after those seen in the 1MB and Kamioka 
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detectors 



• The 5 Baksan events were spread evenly over the entire 10 s period, not bunched 
strongly near the start of the 1987A burst as with the other detectors. Thus, there is 
considerable uncertainty over the true t = time for these events relative to the other 
two detectors. 

The basic idea of the analysis is straightforward. A neutrino of mass m and energy E 
that is emitted at a time tgm and observed in a detector at time t will take a time to reach 
Earth given by 

At = t-tem^toil + m'^/2E^) (1) 

where the approximation requires that |m^| be small compared to E"^ and where to is the 
light travel time from SN 1987A or 168,000 years. Eq. 1 allows us to calculate a mass 
for the kth neutrino event observed at a time tk whose energy is Ek, assuming the neutrino 
was emitted at time tem,k, and travelled a distance cto,fc to reach Earth. Note that to,fc can 
vary from one event to another depending on the collapsing core radius at the time that 
neutrino was emitted. On solving for the mass attributed to a given event, we have: 

2 2E'^{tk — tem,k ~ ^0,fc) 

^k = 7 (2) 

'^0,k 

but clearly eq. 2 can only be used to find rrik when tem,fc and to,fc are both nearly event- 
independent, which requires near-simultaneous emissions. 



II. EVIDENCE FOR NEAR-SIMULTANEOUS EMISSIONS 

The evidence for near simultaneous emissions comes from supernova core collapse models 
and their calculated neutrino fluxes as a function of time, specifically: 

1. brief "burst" of z/g The literature on electron neutrinos and antineutrinos emitted 
during a core collapse consistently shows that both fluxes rise and fall over a very 
short time interval. Typically, this initial burst is found to rise and fall by almost 
an order of magnitude over a time interval of 0.3 seconds, pj, 18[ while some models 
show this "burst" of electron neutrinos and antineutrinos being as short as about 0.02 
seconds, jol 
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2. greater v^. than flux during burst These core collapse simulations find that the 
drop in emitted muon and tau neutrinos following its peak is much more gradual and 
extends over 10 — 15sec, but the initial z/g flux during its short "burst" is found to be 
about an order of magnitude greater than the flux for t > Is,]^ and in any case, 
the SN 1987A neutrinos observed are unlikely to be these flavors - certainly not for 
the 1MB detector whose threshold was below muon neutrino detection. 

3. possibility of brief burst for Although some core collapse simulations find that 
the drop in emitted muon and tau neutrino fluxes following its peak is much more 
gradual than z/g, other authors suggest that a similar rapid rise and fall applies to 
these two flavors as well. lOj - see Keil, et.al for a summary of predictions from various 



models. 
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4. correlation between neutrino energy and emission time If the later arriving neutrinos 
are simply late because they were emitted at later times, one would expect them to 
have higher energies than the earlier emitted neutrinos according to the models, but 
this is not what the data show. 

Given the above, it is not unreasonable to assume that most (or all) of the SN1987A 
neutrinos detected were emitted within & /St ^ 0.3s time interval. Moreover, given the 
rapidity of the core collapse (whose radius drops by an order of magnitude or more depending 
on mass) within a time of 0.2 sec,[l^ we may assume that the neutrino travel distance 
changes very little for those emitted at the start and end of the 0.3s burst. Hence, taking 
tern = ^ ^0) may plausibly assume that within roughly a ±0.2 sec uncertainty tem,k + 'to,k ~ 
iem + to = ± 0.2sec. 

Thus, within this time uncertainty which we assign to the observed t^s, Eq. 2 becomes: 

ml = ^ (3) 
''to ^ ^ 

We note that according to Eq. 3 t = for any E requires m = 0. However, for a pair 
of unsynchronized detectors, for which t = was arbitrarily chosen to mark the arrival of 
their first event, it is unlikely that this time corresponds to the instant a photon emitted 
with the neutrino burst would have reached the detector - especially when the number of 
recorded events in each detector is small. However, straight line fits of versus tk show 
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FIG. 1: On a plot of E~'^ versus t, the SN 1987A data all lie close to one of two straight lines 
that nearly pass through the origin. Solid circles show Kamioka events, triangles 1MB events, and 
open circles Baksan events, which were not used in the analysis to find masses indicative of two 
distinct masses. 



that the data is consistent with one of two lines both cosnsistent with passing through the 
origin. Therefore, apart from a small time offset for each detector t = is consistent with 
tem = ~to- For our estimate of the size of this offset in each detector we have chosen 1/4 
the time difference between each detector's first and second events: 0.107 s in Kamioka, and 
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0.42 s in 1MB in using Eq. 3 to find the mass for eacli event. However, since tliese estimated 
detector offsets are rather uncertain, the ±0.2s uncertainty due to a spread in emission 
times at the source has been increased to a net uncertainty for all the of at = 0.5s when 
computing the uncertainty in the mass associated with each event using: 

^(ml) = ^^\El) + a\t,) (4) 

and CTfc = a{m\) /2mk 



Event 


E in MeV 


t in s 


mass 


K7 


36.9 ±8.0 


1.541 


28.1 ± 5.6 


K8 


22.4 ±4.2 


1.728 


18.1 ± 3.2 


K9 


21.2 ± 3.2 


1.915 


18.0 ± 2.6 


KIO 


10.0 ±2.7 


9.219 


18.7 ±4.5 


Kll 


14.4 ± 2.6 


10.43 


28.6 ±4.8 


K12 


10.3 ± 1.9 


12.44 


22.3 ±3.8 


12 


37.0 ±9.25 


0.42 


14.7 ±3.6 


13 


40.0 ± 10.0 


0.65 


19.8 ±4.6 


14 


35.0 ±8.75 


1.15 


23.1 ±5.3 


15 


29.0 ± 7.25 


1.57 


22.3 ±5.1 


16 


37.0 ±9.25 


2.69 


37.3 ± 8.4 


17 


20.0 ±5.0 


5.01 


27.5 ±6.2 


18 


22.4 ±4.2 


1.728 


18.1 ±3.2 



= 17.8 average mass = 21.4 ± 1.2 



TABLE I: Masses in eV/c^ for 6 Kamioka events and 7 1MB events that are clearly clustered about 
a weighted average of 21.4 eV/(? . Each event's mass is calculated from Eq. 3 using their observed 
time and energy with their associated propagated errors. 

Tables I and II show the results for the 20 events recorded in the two detectors including 
12 H ,n Kanaka. a.d 8 in IMB^ Fo. dar.ty. we have separated the 20 events based on 
whether they correspond to a mass above (Table I) or below (Table II) lOeV^/c^. There 
is a very clear clustering of events into two groups with 13 of them all clustering about 
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mn — 21.4 ± 1.2eV/c^ and the 7 remaining ones clustering about rriL — 4.0 ± O.SeV/c^. 
There is no ambiguity in terms of the cluster to which any one event belongs, since in 
no case was any event less than 2.9(j away from the center of the "wrong" cluster, and 
a majority of events were more than 4(7 away. Moreover, with each cluster the spread in 
event masses was very consistent with what one would expect given the size of the mass 
uncertainties. Thus, the chi squares for all the events fitting one of two masses were quite 
acceptable: xl = 5.9 (6 dof. p=0.43) and x]i = 17-8 (12 dof, p=0.12). In contrast, if the 
20 events were assumed to be described by a single mass, the best value would be 7.9 ± 0.5, 
with — 203.4 (19 dof, p < 10~^°). One might object to this comparison on the grounds 
that an alternative null hypothesis to clustering about two masses would be simply a flat 
distribution of masses, which would in fact give an acceptable chi square. However, given 
that it is known that at most three active mass eigenstates are present, such a distribution 
is physically unrealistic, assuming that the supernova core collapse models are correct in 
predicting a very small spread in emission times. 

The ideogram of Fig. 2 offers a visual way to display the degree of clustering of the 20 
events about two specific masses. In this display each event is represented by a normalized 
Gaussian G{mk, (Jk), with mean and standard deviation and cr^. In addition to showing 
the masses of each event with error bars, we also display the value of ^l^iGlruk, Ck)- 



Event 


E in MeV 


t in s 


mass 


11 


38.0 ±9.5 





0.0 ± 7.3 


Kl 


21.3 ±2.9 





0.0 ±4.7 


K2 


14.8 ± 3.2 


0.107 


3.0 ±1.7 


K3 


8.9 ±2.0 


0.300 


3.0 ±1.0 


K4 


10.6 ±2.7 


0.324 


3.7 ± 1.2 


K5 


14.4 ± 2.9 


0.507 


6.3 ± 1.4 


K6 


7.6 ±1.7 


0.686 


3.9 ± 0.9 



= 5.9 average mass = 4.0 ± 0.5 



TABLE II: Masses in eV/c^ for 6 Kamioka events and 1 1MB event that are clustered about a 
weighted average of 4.0 eV/c^. Each event's mass is calculated from Eq. 3 using their observed 
time and energy with their associated propagated errors. 
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III. DISCUSSION OF RESULTS 



In terms of the robustness of these results, it should be noted that: 

1. Every Kamioka or 1MB event in the data is included in this analysis, using the pub- 
lished values for t, E and Ai?,!!^, Q. Thus, Baksan aside, it would appear there is 
no noise or background events in the SN 1987A data. Moreover unlike earlier inter- 
pretations of the data where only events with t < 2 — 3 seconds are assumed to be 
part of the main supernova burst here it is assumed that they all are. 

2. Every event is cleanly associated with only one of the two mass eigenstates given its 
error bars, and as noted previously the chi squares for the two fits are quite accept- 
able. In contrast fitting all 20 events using a single Gaussian distribution gives an 
infinitessimally small probability. 

3. Every event is consistent with values that are positive. In any case, no assumption 
was made at the outset in this analysis about whether there were or were not events 
having "tachyonic" or superluminal masses, which would have to lie on negatively 
sloped lines in Fig.l. 

4. As a general rule it is not very unlikely to find some kind of interesting pattern with 
only 20 events just based on chance. However, finding a highly specific pattern (every 
single event strongly clustering about one of two mass values is quite another matter, 
especially given the support from supernova models showing that the detection times 
are very likely indicators of the neutrino travel times within a ±0.5s uncertainty in tk- 
One could choose to disregard the present models of supernova core collapse, and insist 
on a significant spread in neutrino emission times. But in that case, one would need to 
assume the existence of some kind of weird correlation between the varying emission 
times and the emitted neutrino energies that perfectly mimics the presence of two 
discrete mass eigenstates at Earth's particular distance from SN 1987A. Moreover, this 
correlation would fail to produce that result if the distance, Lq to SN 1987A happened 
to be significantly different, because the neutrino travel times for a given energy are 
proportional to Lq, making this alternative explanation extremely improbable. 
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FIG. 2: Ideogram of the 20 SN 1987A events computed from the sum of 20 normahzed Gaussians, 
i.e., T,fLiG{mk, (Tk) which displays the clustering of the events around two distinct masses. 



IV. IMPLAUSIBLE MASSES, A MISSING THIRD EIGENSTATE AND OPERA 

The values reported here for the two observed mass eigenstate might seem implausibly 
high. Given that the electron neutrino presumably consists of a mixture of the light mass 
eigenstate plus one or both of the other two, the mass seemingly would need to exceed that 
of ttll = 4.0±0.5e\^/c^, which is inconsistent with the upper limit on the mass of the electron 



neutrino from tritium decay: m < 2eV/(?l^. However, there is an interesting possible 
explanation of this anoma 
neutrinos being correct. 



y_that depends on the OPERA experiment on superluminal ra < 



■ 13] Given the reported OPERA excess over light speed, one can 
calculate the muon neutrino mass that would give that result. Using the average neutrino 
energy in the experiment E = 17GeV one find a superluminal neutrino masss: 

= - v^/c^) = -14, SOOMel/Vc^ (5) 

Since OPERA only measures a single average mass even if several mass eigenstates are 
present, this need not equal the mass of one single eigenstate. However if this mass did 
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represent that of a single superluminal very heavy eigenstate, a SN 1987A neutrino having 
this mass and an energy E ^ 20MeV it would have reached Earth in the past by: 



c 1^1 — m^/E'^ 



which could explain the absence of this eigenstate in the SN 1987 data, i.e., why the 
OPERA result need not be in conflict with the more stringent upper limit from SN 1987A 
on the value of {v — c)/ c. 

Furthermore, the OPERA result, if correct, offers a way to understand why the upper 
limit on the mass of the electron neutrino from tritium decay: m < 2eV^/c^ need not be 
in conflict with the "implausibly high" values of mfj and mi reported here. For example, 
suppose that the electron neutrino consisted of a superposition of a tachyonic and tardyonic 
mass states: \ip >= sin6\%l)^ > + co?,6\il)^ > where the superscripts refer to the signs of the 
respective m? values for the two states whose masses are m? = m\ > and = — < 
To find the mass of the mixed state lip > we need the expectation value of = E"^ — p"^, 
i.e., 



m 



2 



sin9 < tp^\ +COS0 < il)~\{E'^ -p^)sin6'|^+ > +cos6'|V'" > (7) 



If we assume that the mass states have a common momentum p their respective energies 
are E"^ = p^ + and E"^ = — where we obviously have E"^ > p^ > E"^. Thus, Eq. 
7 gives a mass for the mixed state: 

= sin^OEl + cos^9E^__ - p'^ (8) 

It is clear that could be anything between the superluminal (tachyonic) and 
the tardyonic depending on the values of the mixing angle 6. In particular, under this 
scenario, the large superluminal mass eigenstate implied by OPERA could be consistent 
with a value for the electron neutrino mass very close to zero and superluminal depending 
on the mixing angle. In fact, there is previous work that supports the idea that m^g ^ 



-o.25^evyc^.[mm\ 



In the same manner, the possibility of a very heavy superluminal mass eigenstate need 
not conflict with the tight limits from cosmology and astrophysics on the sum of the three 
neutrino masses assuming they have > 0, [21], nor would it necessarily conflict with a 



10 



suggested tachyonic mass limit of that applies to the flavor eigenstates.[22] 

Given the much greater sensitivity of today's neutrino detectors compared to those of 1987, 
were there to be a once-in-a-century occurrence of a supernova in our galaxy, it could offer 
either a clear confirmation or refutation of the existence of the two mass eigenstates claimed 
here. Furthermore, in the event of a confirmation, it would be exceedingly difficult to see 
how the mass limits from tritium beta decay and cosmology could be satisfied without the 
existence of a superluminal mass eigenstate. 
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